Strain driven onset of non-trivial topological insulating states in Zintl Sr2X 

compounds (X=Pb, Sn) 
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We explore the topological behavior of the binary Zintl phase of the alkaline earth metals based 
compounds Sr2Pb and Sr2Sn using both standard and hybrid density functional theory. It is found 
that Sr2Pb lies on the verge of a topological instability which can be suitably tuned through the 
application of a small uniaxial expansion strain (> 3%). The resulting non-trivial topologically 
insulating state display well-defined metallic states in the Sr2Pb(010) surface, whose evolution is 
studied as a function of the film thickness. 

PACS numbers: 71.20.-b, 73.43.-f, 73.20.-r 
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Topological insulators [l|-|6[ (TI) are a new quantum 
state of matter characterized by the existence of gap- 
less surface states sheltered against destructive scattering 
effects by time-reversal symmetry. Since the discovery 
of two dimensional TI behaviors of HgTe-based quan- 
tum wells several families of topological materi- 
als have been theoreti cally predicted and experimentally 
realized afterwards. [ll The variety of insulating ma- 
terials displaying topological features includes BixSbi^x 
allovsjR Bi2Te3, Bi2Se3 and Sb2Te3 binary com- 
pounds El \lM > ternary heayy-metals based com- 
20unds such as TlBiTea, TlBiSe2 ^^E^, and PbBi2Se4 
17,(3], and ternary rare-earth chalcogenides (LaBiTca) 
19| , and another honeycomb-lattice type of ternary com- 
pounds (LiAuSe) [2^. All these materials are character- 
ized by a layered structure stacked along the c-axis of the 
centrosymmetric hexagonal lattice, similar to the struc- 
ture of Bi2Se3. In addition, other ternary TIs have been 
recently predicted such as the non-centrosymmetic cu- 
bic zinc-blende HgTe-like phase (half-Heusler compounds 
[HKii and I-III-VI2 and II-IV-V2 chalcopyrite semicon- 
ductors such as AuTlTe2 [3, IH), and CasNBi with a 
centrosymmetric antiperovskite structure |26[ . Although 
a large number of ternary TI have been found 3- 3 21 - 
[2^, up to now, binary TIs are only limited to three 
classes: Bii_2;Sba; alloys and the family of Bi2Te3, Bi2Se3 
and Sb2Te3 as well as Ag2Te [27|]. In searching for new TI 
materials we have focused our attention on binary heavy- 
element-based small band gap semiconductors. One of 
the simplest way to realize a binary heavy-element-based 
closed-shell semiconductor is M2X with an alkaline earth 
element (M ~ Mg, Ca, Sr, and Ba) which donates its s 
valence electron to a IV group element {X ~ Si, Ge, Sn 
and Pb). These are the so-called Zintl compounds. 
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FIG. 1: (Color online) Structure representations of or- 
thorhombic (space group Pnma, No. 62) M2X compounds 
(M = Mg, Ca, Sr and Ba; X = Si, Ge, Sn and Pb): (a) the 
unit cell (u.c), (b) the projection of the unit cell in the ac 
plane, perpendicular to the 6-axis and (c) the local environ- 
ment around the X atom. The large (gray) and small (red) 
balls denote the M and X atoms, respectively. In addition, 
along the &-axis (y=l/4 and y=3/4) the atoms can be ar- 
ranged in two parallel planes shown in the dashed and solid 
balls, respectively. M occupied two inequivalent 4c sites, M\ 
[xi, 1/4, zi) and M2 {x2, 1/4, Z2) whereas X occupies the 4c 
site (xa, 1/4, 23). For Sr2Pb the optimized atomic positions 
are Sri (0.0214, 1/4, 0.683), Sr2 (0.1580, 1/4, 0.073) and Pb 
(0.2496, 1/4, 0.3929). 



The Zintl compounds M2X crystallize in a simple or- 
thorhombic crystal structure [13, [2^, as illustrated in 
Fig. [1] The X atoms surround the M atoms in a slightly 
distorted trigonal-prismatic coordination. The coordi- 
nation sphere is augmented by the three X atoms situ- 
ated on the rectangular faces of these prisms (tri-capped 
trigonal prism). Along the a direction the prisms are 
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FIG. 2: (Color online) DFT and HSE electronic structure and band inversion strength. Panels (a,b) and (c,d) show DFT 
calculated electronic band structures for Sr2Sn and Sr2Pb, respectively. Panels (e,f) and (g,h) reports HSE band structures 
around the F point along the Z-F-T directions for Sr2Pb at zero strain and under 5% tensile strain in the ac-plane, respectively. 
The upper and lower panels show the results calculated without (a,c,e,g) and with (b,d,f,h) SOC, respectively. The solid (red) 
circles denote the states with a predominant s-like character. Panel (i) illustrates the comparison of the band inversion strength 
between p^- and s-orbitals at the F point as a function of strain. Negative values denote the occurrence of band inversion. 



condensed via common edges while they share common 
triangular faces along the b direction. The semiconduct- 
ing character of these compounds originates from their 
closed-shell nature, with 6 valence electrons (2xs^+p^) 
per formula unit [{2]\/P+)X'^-]. 

We have performed band structure calculations using 
the Perdew-Burkc-ErnzcrhoffsO] (PBE) based standard 
and hybrid (HSE^]) density functional theory (DFT) as 
implemented in the Vienna Ah initio Simulation Package 
(VASP) [13, IISIj with the inclusion of relativistic spin- 
orbit coupling (SOC) effects. Within the HSE method, 
the many body exchange and correlation functional is 
constructed by mixing 25% exact Hartree-Fock exchange 
with 75% PBE, and the long-range Coulomb interaction 
is suitably screened according to the parameter ^ (here /i 
= 0.3 A-i). By using the experimental lattice constants 
we have relaxed all atomic positions using a very dense 
k-points mesh (up to 4200 kpoints). 

Our preliminary PBE-based research of possible topo- 
logical features in this class of materials revealed that two 
members of the Zintl family (Sr2Pb and Ba2Pb) exhibit 
promising fingerprints of TI behaviors whereas the re- 
maining compounds are either trivial insulator (Ca2Ge, 
Sr2Ge,Ca2Sn and Sr2Sn) or semimetal (Ba2Ge, Ca2Pb 
and Ba2Sn). Therefore in the following we focus our 
analysis on these two specific Zintl compounds: Sr2Pb 
and Sr2Sn. 

The DFT band structures with and without SOC ef- 
fects are compared in Fig. [5l^a-d). The inclusion of 



SOC effects does not affect the overall electronic char- 
acter of Sr2Sn [Fig. [5] (a,b)] which remains semiconduc- 
tor with a small bandgap at F (« 0.15 eV) opened be- 
tween occupied Sn states and highly dispersive Sr d- 
like empty orbitals (mostly d^2_j^2 ). By replacing Sn with 
the isoelectronic heavier element Pb the band structure 
changes dramatically, as schematically depicted in Figl3] 
(a). Without SOC, Sr2Pb displays a metaUic (gapless) 
state [Fig. [2][c)] whereas the inclusion of SOC [Fig. [UJd)] 
opens a band gap of about 100 meV at F (the indirect gap 
is about 50 meV), as a consequence of the anti-crossing 
between the conduction band minimum (CBM) and va- 
lence band maximum ( VBM) , a typical fingerprint of the 
SOC-induced formation of topological insulating states. 
The major differences between the electronic structure of 
Sr2Sn and Sr2Pb resides in the orbital character of the 
valence and conduction bands near F as highlighted by 
the (red) solid circles in the band plot of Figl^Ja-d): in 
Sr2Sn the states with a predominant s-like character lie 
about 0.5 eV above the Fermi Level [see Fig. [2][a,b)], 
whereas in Sr2Pb these states are pushed down in energy 
and eventually hybridize the VBM at F thus inducing 
the anti-crossing feature responsible for the opening of 
the gap and the creation of a TI state. The downward 
shift of the s-states is accompanied by a upward shift of 
the heavy metal p bands which ultimately intermix with 
the Sr d states as schematized in Fig. 131(a). 

As recently reported by Zunger and coworkers [3^. the 
identification of band-inverted Tl[^, [l3| on the basis of 
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conventional DFT may lead to false-positive assignment 
if the band inversion strength is not large enough, due to 
the well documented bandgap underestimation problem. 
For this reason we have revisited the electronic dispersion 
of Sr2Pb by HSE, which indeed yields a substantially dif- 
ferent physical picture, as illustrated in Fig. [Ifejf). HSE 
finds a much larger bandgap (0.25 eV) which is reduced 
to 0.13 eV with the inclusion of SOC effects, and, most 
importantly, prevent the occurrence of band inversion be- 
tween the pa;-like and s-like states as schematically illus- 
trated in Fig. ^h). The s states remain well localized on 
the bottom of the conduction band [Fig. [2je,f)]: Sr2Pb 
is not a TI in its native phase. DFT wrongly stabilizes a 
spurious TI solution because of the relatively small band 
inversion strength [0.12 eV, see FiglSfa)]. Though HSE is 
expected to provide a generally more accurate description 
of band dispersion in small bandgap insulators, future op- 
tical experiment will be necessary in order to validate our 
first principles findings. 
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FIG. 3: (Color online) Evolution of atomic orbitals at the F 
point from Sr2Sn to Sr2Pb with and without SOC effects in- 
cluded for both DFT (a) and HSE (b) calculations. The band 
inversion can be observed by Pb (or Sr) s and Pb Px orbitals. 
In panels (c) and (d) we show the product of wave function 
parities of the occupied bands for eight time-reversal invariant 
momenta (TRIM) in the Brillouin zone [F (0,0,0), X (7r,0,0), 
Y (0,7r,0), Z (0,0,7r), S (7r,7r,0), T (0,7r,7r), U (7r,0,7r), and R 
(7r,7r,7r)] of Sr2Pb obtained before and after band inversion, 
respectively. 

To explore possible routes for designing a TI phase 
in the native phase of Sr2Pb, we carried out a series 
of calculations applying strain, a gap-engineering tech- 
nique which was successful when applied to zero-gap 



semiconductors such as ternary Heusler compounds [2^ 
and CaaNBi. [2^ We have chosen to study the effect of 
uniaxial strain (e) in the ac-plane, by leaving the 6-axis 
unconstrained (free to relax) in order to simulate at best 
the experimental condition for thin film growth. The re- 
sults, obtained by HSE and shown in Fig. [SJg and h), 
indicate that for relatively small uniaxial strain larger 
than 3% Sr2Pb can be tuned towards a TI phase. The 
role of SOC effects is essential to open a small gap around 
F and to induce an inverted band order. The s-like states 
shift downward below the Fermi level and becomes occu- 
pied and, simultaneously, the Pa:-like states become un- 
occupied and promoted at higher energy. This kind of 
inverted band behavior can be ascribed to the fact that 
the strain-induced expansion in the ac-plane reduces the 
crystal field effect, resulting in the less hybridization be- 
tween s- and p-like states and stabilizing the s-like state 
at lower energy, as evidenced in Fig. ^e) and (g). In such 
a way, the spin-orbit coupling strength is now enough to 
invert the band order between s-like and p-like states at 
r (c.f., Fig. [2Jf) and (h)). For larger strains up to 7% the 
TI state is preserved and further stabilized as inferred by 
the evolution of the band inversion strength as a func- 
tion of strain reported in Fig. [5Ji) at both, DFT and 
HSE levels. 

Beside the band inversion, an alternative way to iden- 
tify TI states is the parity criteria proposed by Fu and 
Kane I4I . Considering that orthorhombic Pnma posses 
the inversion symmetry this criteria can be applied and 
will serve as a further support for our analysis. The prod- 
uct of the parities of the Bloch wave function for the 
occupied bands at all eight time reversal invariant mo- 
menta (TRIM), illustrated in Fig. [2Ic,d), suggest that 
at six TRIMS (X, Y, Z, S, T and U) ah bands share the 
same doubly degenerate character, whereas the TRIM R 
is found to be fourfold degenerate. Therefore all these 
seven TRIMs display a positive ("-I-") global parity. At 
r the situation is different: the product of the parities is 
and "-" depending on whether or not the band in- 
version occurs, in consistency with the HSE band struc- 
ture interpretation [Fig. [3l[c,d)]. We can therefore trust- 
fully conclude that distorted orthorhombic Sr2Pb is a 
topological non-trivial insulator with Z2 index (1; 000). 
Strained-driven gap-engineering on Sr2Sn dose not re- 
sults in any topological transition: Sr2Sn remains a con- 
ventional semiconductor. 

After discussing the onset of topological features in 
the bulk phase of distorted Sr2Pb we turn our atten- 
tion to the surface properties focusing on the non-polar 
(010) termination. Considering that both DFT and HSE 
lead to an essentially identical TI state in strained Sr2Pb 
and that SOC-HSE calculations for thick slabs are com- 
putationally very demanding (if not prohibitive at all) 
we study the surface band structure at DFT level only. 
The results on the Sr2Pb(010) surface are summarized 
in Fig. 21 which shows that very thick slabs have topo- 
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logically protected surface metallic states which remain 
robust with increasing the film thickness, thus corrobo- 
rating the conclusion of the bulk parity analysis and the 
inverted band order. The evolution of the band structure 
as a function of the film thickness show a very peculiar 
behavior. Already at low film thickness surface related 
states emerge in a small energy window (± 200 meV) 
around the Fermi level. This leads to a quenching of 
the bandgap as compared to the bulk value (130 meV). 
The bandgap is then progressively reduced by increasing 
the slab thickness. Though for a thickness of 22 unit cell 
(u.c.) [Fig. mjf)] the gap is almost completely suppressed 
(10 meV), it is necessary to increase the thickness up to 
100 u.c. in order to fabricate a well-defined metallic film 
[Fig. m^g)]. At the critical thickness of 100 u.c. the low- 
est conduction band and the highest valence band meet 
at the Fermi level at F and establishes a single-Dirac- 
cone-like metallic surface state. 
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FIG. 4: (Color online) DFT Surface properties of strained 
Sr2Pb(010) (e =5%). (a-d) Evolution of the band structures 
with slab thickness (e) structural model of the symmetric slab 
adopted to simulate the Sr2Pb(010) surface (the image cor- 
respond to the 3 unit cells case, i.e. 6 layers per side. We 
adopted a vacuum region of 30 A), (f) Band structure for a 
thickness of 22 u.c. Shaded areas refer to bulk bands, (g) Pro- 
gressive closing of the gap at F (E^) as a function of the slab 
thickness. Note that, due to the prohibitive computational 
cost, for slab thickness larger than 34 u.c. (> 400 atoms) the 
band gaps were calculated at F only. 

In conclusion, our computational study has disclosed 
the non-trivial topological nature of the binary com- 
pound Sr2Pb. The detailed analysis of the bulk and sur- 
face structural and electronic properties is of relevance 
for the design principles of TI and provides helpful in- 
sight for tunability of TI states in trivial insulator by 



gap-engineering techniques. We believe that our findings 
will encourage immediate experimental investigations. 
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